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ABSTRACT 


Previously published methods of measuring and interpreting the Hall coefficient and 
resistivity of semiconductors are summarized. The conventional and van der Pauw 
methods of measurement are briefly discussed. Mobility is shown to be a function of 
carrier scattering in a specific semiconductor. The equations for charge carrier con- 
centration are developed for both extrinsic and intrinsic conduction. These equations 
are then used as the basis for interpreting the measured experimental data. The analy- 
sis of the data near the intrinsic range and the temperature range over which the experi- 
mental measurements must be made are specified. 
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A SUMMARY OF THE MEASUREMENT AND INTERPRETATION OF THE 
HALL COEFFICIENT AND RESISTIVITY OF SEMICONDUCTORS 
by D. Michael Stretchberry 
Lewis Research Center 

SUMMARY 

A knowledge of the parameters of semiconductor material is important both for basic 
materials research and for device development. An analysis of Hall coefficient and re- 
sistivity data measured over a temperature range will yield several of these parameters, 
such as the impurity doping concentrations, the impurity ionization energy, the band gap 
energy, and mobility. The measurement and interpretation of the Hall coefficient and re- 
sistivity are summarized from the available literature. 

The conventional and van der Pauw methods of measuring Hall coefficient and resis- 
tivity are briefly discussed. The van der Pauw method was found to be advantageous for 
small samples. 

The mobility of semiconductors is also discussed. It is related to an energy- 
averaged relaxation time, which is a characteristic of a particular scattering process and 
a specific semiconductor material. The different kinds of scattering are qualitatively 
discussed. The factor a, which relates the Hall coefficient to the carrier concentra- 
tion, is also a function of relaxation time. 

The equations for charge carrier concentration are developed for both extrinsic and 
intrinsic conduction. They are found to be functions of temperature and the various semi- 
conductor parameters. These equations are then used as the basis for interpreting the 
measured experimental data. The method of interpretation is to estimate the equation 
parameters and then to calculate the carrier concentration at each temperature of inter- 
est. The parameters are varied until the calculated results match the experimental data. 

The shape of the Hall coefficient and resistivity against temperature curves near in- 
trinsic conduction is shown so that intrinsic data can be interpreted correctly. 

The temperature range over which the experimental measurements must be made is 
discussed. There are certain characteristic parts of the plot of the logarithm of the car- 
rier concentration against reciprocal temperature from which data must, be included for 
extrinsic or intrinsic interpretation. 



INTRODUCTION 


The NASA- Lewis Research Center is currently engaged in a program which includes 
the determination of the basic properties of semiconductor materials, particularly for 
high-temperature applications. The ultimate objective is the development of devices that 
operate at temperatures greater than 500° C (773 K). The measurement and analysis of 
the Hall coefficient and resistivity as a function of temperature was picked as the method 
for determining, the desired properties. Although there are other ways of obtaining these 
properties, this method is the most popular and has been used since the early days of 
semiconductor research. The derivations of the pertinent equations and the methods of 
analysis which result are well known (refs. 1 and 2). But, it was realized that the liter- 
ature does not emphasize the use of the derived equations for the interpretation of exper- 
imental data and that no single source contains all of the required information. There- 
fore, a literature search was started, and this summary report resulted. The purpose 
of this report is to make available the collected information to the researcher who wants 
to begin using the Hall coefficient and resistivity methods. More details of the methods 
can be found in the references cited. This report has proved very useful to our pi'ogram 
and should be equally useful to others. 

This report summarizes the use of the Hall coefficient and resistivity data as a func- 
tion of temperature to obtain several of the electrical properties of semiconductors, such 
as mobility, doping concentrations, the ionization energy of the major dopant, and the band 
gap energy. The topics discussed are applicable to all semiconductors and not limited to 
those suitable for high-temperature use. The introductory sections are concerned with the 
applicability of Hall coefficient and resistivity data (ref. 3), the measurement of these quan- 
tities (refs. 1 and 4), and the charge carrier mobility in semiconductors (refs. 5 to 8). Thes 
sections are included for a complete discussion of the problem of interpreting the data. 

The applicability of the Hall coefficient and resistivity to a basic materials research 
program and to a device development program is illustrated by specific examples. Sev- 
eral of the parameters upon which device design depends are directly obtained from Hall 
data. The conventional and the van der Pauw methods of measurement are briefly de- 
scribed (refs. 1 and 4). The latter method is valuable because it is independent of the 
shape of the sample and, therefore, is more easily applicable to small samples. The 
mobility of semiconductors is discussed both as it illustrates the scattering processes 
taking place in the semiconductor and in its role of influencing the relation between the 
Hall coefficient and the charge carrier concentration (refs. 5 to 8). The bulk of the dis- 
cussion is qualitative because of the complexity of mobility theory. 

The equations for the charge carrier concentration were developed both for extrinsic 
and intrinsic conduction (ref. 2). For the extrinsic case they are a function of tempera- 
ture, impurity concentrations, and the impurity ionization energy. The oversimplified 
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uncompensated case is considered first and then generalized to the more realistic com- 
pensated case. The charge carrier concentration equation for intrinsic conduction is 
shown to be a function of the temperature and the band gap energy. 

The bases for interpreting the measured experimental data are these equations. The 
carrier concentration is calculated at the required temperatures by using various sets of 
estimated values of the constants. By comparing the calculated results with the experi- 
mental data, those values of the constants that allow the best curve fit can be obtained. 
Initially, the constants of the equations are estimated from the experimental data curves. 
The characteristic shapes of the Hall coefficient and resistivity against temperature 
curves in the temperature range near intrinsic conductivity is illustrated (ref. 9). Knowl- 
edge of these characteristic shapes prevent errors being made in the interpretation of the 
data in this region. 

Concluding the report is a brief discussion of the relevant criteria for specifying the 
temperature range over which the Hall coefficient and resistivity must be measured. 

Data must be obtained over a temperature range wide enough so that the methods of in- 
terpretation that are to be developed can be used. 


APPLICABILITY OF HALL COEFFICIENT AND RESISTIVITY 


A knowledge of the Hall coefficient and resistivity of a semiconductor is important 
both for basic materials research and device design. These measurements have proved 
useful since the discovery of semiconductors, and their value is illustrated by specific 
examples. 

The band gap energy and energy levels of the impurities are basic electrical param- 
eters of the material and are valuable in a semiconductor materials research. Also 
valuable is information about the transport properties of the charge carriers which an 
analysis of the mobility data against temperature gives. 

Semiconductor device design depends on several parameters, such as the impurity 
concentrations (doping levels) of each material, the physical dimensions of each part of 
the device, the basic constants particular to the device material, the diffusion constants 
(or mobilities) of the charge carriers, the lifetime of the carriers, and the proposed op- 
erating temperature of the device. A knowledge of the Hall coefficient and resistivity as 
a function of temperature gives the designer the doping levels of the materials, the car- 
rier mobilities (and therefore the diffusion constants), and the band gap energy. The dif- 
fusion constants can be obtained directly from the mobilities by Einstein's relation. The 
band gap energy affects the operating temperature of the device. The basic constants of 
the material (such as the breakdown electric field or the dielectric constant) can usually 
be measured. Normally, the physical dimensions are measured or calculated. The life-. 
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time of the charge carriers is measured by methods not related to the Hall coefficient or 
resistivity. Therefore, it is evident that a substantial amount of information (but not all) 
is ihade available by an analysis of the Hall coefficient and resistivity. 

The single- junction (diode) design equations are dependent, in part, on the doping 
levels and diffusion constants. For example, the junction capacitance of an abrupt junc- 
tion is proportional to the square root of the dopant density of the more lightly doped re- 
gion. The breakdown voltage is inversely proportional to this dopant density. The as- 
sumption is made that the doping atoms are completely ionized. The diffusion constant 
appears in equations pertaining to current flow. Two examples are the characteristic 
voltage- current equation of the diode and the expression for the ac impedance. 

Other devices are made by combining junctions in certain configurations. There are 
specific considerations for each type of device, but a common requirement is good junc- 
tion design. For transistors, the current gain and the four pole parameters (which spe- 
cify the ac characteristics) are both functions of doping levels and diffusion constants. 

It is evident that both the basic materials researcher and the device developer are 
interested in measuring and interpreting Hall coefficient and resistivity data. The in- 
formation obtained from these relatively simple electrical measurements is of much 
value. 


MEASUREMENT OF HALL COEFFICIENT AND RESISTIVITY 


A brief description is given of both the conventional and van der Pauw methods of 
measuring the Hall coefficient and resistivity. References are cited to indicate where 
comprehensive information about both methods can be found. 


Conventional Method 

The conventional Hall and resistivity measurements are made on constant-thickness, 

single- crystal samples of semiconductor material, as shown in figure 1. The arms on 

/ 

each Side are for ease in making contacts which will perturb the current flow through the 
length of the sample as little as possible, A current I (in amperes) is passed the length 
of the sample from A to B. If the voltage drop (in volts) from C to E is measured 
and the dimensions of the sample (in centimeters) are known, the resistivity can be cal- 
culated as follows: ’ v • - 


v CE tw 
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Figure 1. - Bridge-shaped sample for conventional hall and re- 
sistivity measurements. (B peroendicular to plane of crystal.) 


where t, W, and L are particular dimensions of the sample as shown in figure 1. The 
conductivity (which is the inverse of the resistivity) can usually be written as 

n-type (electrons majority carrier): 


a = 


1 

P 


n o^n 


p-type (holes majority carrier): 


( 2 ) 


( 7 


1 

P 


p 0 qp 
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where n Q and p Q are the equilibrium concentrations, ju n and /i^ are the conductivity 
mobilities of electrons and holes, respectively, and q is the electronic charge. If both 
electrons and holes are present in equal or nearly equal concentrations, the conductivity 
equation is the more general equation 


a 


i 

p 


<l(Vn + ¥p> 


(3) 


Analogous equations can also be written for the Hall coefficient. If in addition to the 
current I a known magnetic field B (in tesla) is applied perpendicular to the plane of the 
crystal, a voltage V H will be observed across the width of the crystal from C to D. 

A measurement of V H will allow the calculation of the Hall coefficient R^. It is ex- 
pressed in terms of measurable quantities as 


R H ~ 


_V 


BI 


•xlO 4 cm vC 


(4) 
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The sign of Rjj indicates the material type (i. e. , for R^. > 0, the semiconductor is 
p-type; for Rg < 0, it is n-type). 

The Hall coefficient is 

n-type: 



p-type: 


( 5 ) 



where a is a factor with a value usually between 1 and 2. The factor a is discussed later 
(in the section MOBILITY AND ITS INFLUENCE). If both types of carriers are present in 
equal or nearly equal concentrations, the Hall coefficient is more generally expressed as 



p - n 
*Vp o'n 

(P« Ik* + 

Vi o p o'n 7 


( 6 ) 


It can be seen from this equation that, if the Hall coefficient can be measured, the 
carrier concentration can be easily calculated for the case where there is only one type 
of carrier. When there are two types of carriers, the concentrations are not easily de- 
termined. This is discussed later in the section INTERPRETATION OF DATA. The mo- 
bility can be calculated when the Hall coefficient and resistivity are known, as is discussed 
in the section MOBILITY AND ITS INFLUENCE. 

The material presented in this section can be typically found in elementary books on 
semiconductors such as Shockley’s (ref. 27). 


van der Pauw Method 

A method of measuring the Hall coefficient and resistivity of arbitrarily shaped sam- 
ples has been derived by van der Pauw (ref. 4). Only four contacts instead of the five or 
more normally required are necessary for both measurements. The most important ad- 
vantage of this method is obvious: the sample no longer has to be machined to a particu- 
lar shape. This is advantageous when only small crystals are available and machining 
would waste a large proportion of the sample. The only restrictions of van der Pauw’s 
derivation are (1) the sample must be of constant thickness, (2) there must not be cavities 


6 




Figure 2. - Arbitrarily shaped van der Pauw 
sample with jour electrical contacts (a, b, 
c, andd). (B perpendicular to plane of 
crystal. 


in the crystal, and (3) small contacts must be made as close to the perimeter as possible. 
A van der Pauw sample is shown in figure 2. 

The resistivity equation derived by van der Pauw is 


P 


wt 


2 In 2 


(“ab, cd + “be, da> f 0^') ohm - cm 
Vbc, da/ 


(V 


The function f is given by van der Pauw in graphical form (see fig. 3). The quantity 
R a b cc j is defined, for example, to be the ratio of the voltage measured from d to c to 
the current flowing in a and out b. Other such ratios are defined similarly and are also 
expressed in ohms. 

The Hall coefficient expression defined by van der Pauw is 


R H=j i V,irf !ll)4< ” 3/c < 8> 



Figure 3. - Function used for determining resistivity of sample, plotted as function of 
R ab, cd/ R bc, da- 
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where AR ac ^ is the difference between R ac ^ without the applied magnetic field and 
with the field. The sign of Rjj indicates the type of the semiconductor. 

Influence of Thermomag netic and Thermoelectric Effects 

Thermomagnetic and thermoelectric effects can influence the Hall voltage. Some, 
such as the Seebeck, Nernst, and Righi-Leduc effects, can be eliminated through averag- 
ing the data obtained by reversing the current direction. Another, the magnetoresistance 
effect, is eliminated through averaging the data obtained by reversing the magnetic field 
(refs. 1 and 7). Finally, the Peltier-Nernst and Ettinghausen effects cannot be removed 
by data averaging but usually are small enough to neglect (ref. 7). However, their in- 
fluence should be estimated to verify this assumption. 


Suggested References 

A popular book which describes nearly all aspects of the Hall effect is that by Putley 
(ref. 1). The conventional measurement of Hall coefficient and resistivity, as well as the 
general equipment that is required, are discussed. The thermomagnetic and thermoelec- 
tric effects and their influence on the Hall coefficient are described in detail. 

The original paper by van der Pauw (ref. 4) gives the basic information on his 
method of measurement. In reference 7 van Daal extends the discussion of the van der 
Pauw method and also describes the influence of the thermomagnetic and thermoelectric 
effects. 


MOBILITY AND ITS INFLUENCE 

Several kinds of mobility are used in the semiconductor literature. They are seldom 

t 

differentiated well. Mobility is defined herein since it is used often in this report. When a 
semiconductor is subjected to an electric field, the charge carriers are given a drift ve- 
locity. The drift velocity is proportional to the electric field. The constant of propor- 
tionality is defined as the mobility \i. The basic types of mobility depend on whether the 
carriers are excess minority carriers or majority carriers in equilibrium (ref. 10). 

The mobility for the former case is called drift mobility and for the latter case is called 
conductivity mobility. The distinction between the two is usually not significant. But, 
precisely speaking, the conductivity mobility is the appropriate mobility to use in all 
cases in this paper. The mobility of electrons is generally different than that of holes. 
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It is possible to have electron and hole conductivity mobilities, as well as electron and 
hole drift mobilities. 

There is in general use another quantity called the Hall mobility. It is a type of ‘con- 
ductivity mobility and is defined as 



( 9 ) 


If equations (2) and (5) are used, it can be shown that the Hall mobility is related to the 
conductivity mobility by 


M h = (10) 

Interest in mobility as it relates to this report is twofold. First, the analysis of 
mobility as a function of temperature yields information about the scattering processes 
which influence the movement of charge carriers in the semiconductor. Secondly, a is 
obtained in the study of mobility as a function of the scattering processes and of the par- 
ticular semiconductor band structure. 


Relating Mobility to Scattering 

The relation between mobility and scattering shows that the mobility is directly de- 
pendent on the energy average of the relaxation time r. The energy average of r is de- 
noted as (t) . The relaxation time is closely related to the average time between col- 
lisions of the charge carriers. They are equal only when r is not a function of energy. 
(See Smith (ref. 5) for more details. ) The relaxation time is introduced as the charac- 
teristic time governing the establishment of equilibrium of the initially disturbed carrier 
concentration. The relaxation time as a function of energy (or, more generally, velocity) 
must be calculated for the specific scattering process and then energy- averaged to ob- 
tain the mobility. The factor a also depends on (r) . 

The derivation is based on the following assumptions: 

(1) Spherical constant- energy surfaces ink- space (k- space is related to momentum 

space by a multiplicative constant) 

(2) An n-type semiconductor with a parabolic conduction band 

(3) A nondegenerate semiconductor 

(4) Relaxation time r as a function of energy only (isotropic) 

The first assumption presumes that the properties of the semiconductor (such as conduc- 
tivity and effective mass) are isotropic. The second assumption allows the usual density- 
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of- states expression to be used in the derivation (see the section DEVELOPMENT OF 
CHARGE CARRIER CONCENTRATION EQUATIONS). The third assumption allows the 
use 1 of Boltzmann statistics in place of Fermi-Dirac statistics. Also assumed is an elec- 
tric field in the x-direction £ x . This derivation appears in Smith (ref. 5), Putley 
(ref. 1), vanDaal (ref. 7), and Blatt (ref. 6). 

The result is an equation for J x , the current density in the x-direction, as a function 
of (t) , which yields the results for mobility and the factor a. This current density 
equation is 


J = 


n o q 2 ^ x (r) 

m* 

n 


where m* is the effective mass of electrons and (r) is defined to be 


(r) = 


3? r 


1/2 


/ 

•'O 


exp(-e)de 


( 11 ) 


( 12 ) 


where e is 


e = 


E - E„ 
c 

kT 


(13) 


and k is the Boltzmann constant, T is temperature, E is energy, and E is the energy 
of the bottom of the conduction band. The average drift velocity of the electrons is given 
by 



This relation plus the definition of mobility 




l< T x>l 


X 


(14) 


( 15 ) 


allows the mobility to be expressed in terms of the energy- averaged relaxation time 
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(16) 


M n 


m n 


The definition of conductivity is 



Thus, it can be expressed as 


cr = n Q qp 


n 


(17) 


(18) 


and equation (2) is verified. 

The mobility has now been expressed as a function of the energy-averaged relaxation 
time. To calculate the dependence of mobility on any scattering process all that need be 
done is to determine the characteristic r for that type of scattering and calculate its 
average by use of equation (12). However, this is difficult to do. 

If the same type of calculation as shown in equations (11) to (18) is carried out, but 
including the effect of a magnetic field on the charge carriers, the result for the factor 
a is 


/_ 2 \ 

(19) 

( 7 -) 2 

Generalizing Mobility 

2 

The equations for p, and a. in terms of (r) and (t ) will be modified when the 
constant- energy surfaces are not spherical or when there are several equivalent energy 
minimums in the conduction band. Consider first the case of nonspherical constant- 
energy surfaces with a single energy minimum in the conduction band at k = 0. (The n 
subscript will be dropped from the mobility symbol.) 

The effective mass is a function of the shape of the energy bands and therefore, in 
general, is a tensor quantity. It is scalar only when the constant- energy surfaces are 
spherical. However, when they are not spherical, a set of axes can be chosen to make 
the effective mass tensor diagonal. The applied electric field then can be broken up into 
its components along the chosen axes. 
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By an extension of this derivation, a set of equations similar to equation (11) results: 

"\ 


J x' = 




m; 


_ _ <*>* Y l 
y " mi 


J z = 


n o q <*>-*■ 


m? 


( 20 ) 


J 


where m^, m|, and m* are the diagonal elements of the effective mass tensor, and 
therefore 




H- 

mj 


mS 


ms 


( 21 ) 


The conductivity is also a diagonal tensor in terms of these axes and its components are 
written as 

°i = 


ct 2 = n o q/ *2 ? 


( 22 ) 


= n o qM 3 


Now consider the existence of several equivalent energy minimums in the conduction 
band. If there are M minimums, there are n Q /M electrons per unit volume in each. 
Assume M = 6, which is the case for silicon (see fig. 4). The sum over all the electrons 


12 




Figure 4 - Constant-energy surfaces in k-space for energy Just above 
bottom of conduction band i n sil icon. Mini mums are at centers of el- 
lipsoids. 
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a = n Q qii 




M 


m n 


J 


( 25 ) 


the conductivity effective mass and conductivity mobility are defined as 


m* 

n 



M - (M j * M 2 + M3) 
3 


(26) 


Note that the constant -energy surfaces for silicon are ellipsoids. Therefore, two of 
the effective masses in equations (26) will be equal (m^ = m| = m*,) because of the sym- 
metry of the ellipsoids, and will be different from the third (m* = m£). Equations (26) 
then reduce to 


i = l(_L + JL\ 

m* 3 \ m li m T/ 


IX - — (/X T + 2firp) 

3 u 1 


(27) 


The conductivity effective mass and conductivity mobility are defined because the 
equations in which m^, m|, m*, jx^, p 2 > and M 3 appear can be written more conven- 
iently with m* and fx. It must be remembered, however, that m* and /x depend on 
the band structure of the semiconductor. Thus, equations (26) and (27) are not general 
equations for all band structures. A Hall effective mass (which is analogous to the con- 
ductivity effective mass) can be defined for use with the Hall mobility. 

The factor a is also more complicated for the case when there are several equiv- 
alent energy minimum s in the conduction band and nonspherical constant- energy surfaces. 
If the number of equivalent minimums is six, for example, a is given by (ref. 5) 


a 



3(mJ + m* 



1 m*m*m* 


(28) 
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If the constant- energy surfaces are ellipsoidal, m| = m| = and m| = m£, a be- 
comes 


q - fr 2 ) 3C(C + 2) 
<r> 2 (2C+1) 2 


( 29 ) 


where C is defined as m* /m* 


Scattering Processes 

The derivation of a relaxation time for a particular type of scattering is not presented 
because of the complexity of the calculations. The different types of scattering processes 
are, however, qualitatively discussed. In general, the mobility will be determined at low 
temperatures by crystal structure defects and impurities, and at high temperatures by 
vibrations of the lattice. Table I lists the different possible types of scattering processes 
and the factor a corresponding to them. The values for a are from van Daal (ref. 7). 
They include only the (t )/(t) part of a. The effective mass part is assumed to be 
equal to 1. References are also listed in table I. 

Lattice vibrations are of two types. Acoustic vibrations are those in which an atom 
of a crystal and its nearby neighbors vibrate in the same direction (in phase). Optical 
vibrations are those in which two adjacent particles move in opposite directions (out of 
phase). Transverse and longitudinal vibrations can exist for both types of lattice vibra- 
tion. Figure 5 illustrates acoustic and optical vibration for both the transverse and lon- 
gitudinal cases. 

Acoustic and optical scattering differ for polar and nonpolar crystals (see table I). 
Acoustic scattering for nonpolar crystals is called acoustic-mode scattering and for polar 
crystals is called piezoelectric scattering. Optical scattering for nonpolar crystals is 
called optical-mode scattering, and for polar crystals it is called polar scattering. Sili- 
con and germanium are nonpolar crystals, and silicon carbide is polar. Whether a crys- 
tal is nonpolar or polar is determined by whether the crystal’s bonds are polar or non- 
polar. A polar bond will have the center of negative charge offset from the center of pos- 
itive charge. Ionic bonding is polar, but covalent bonding can be either. 

Early calculations of mobility considered only acoustic-mode scattering. This lead 
to a mobility which varies as T~ 3 / 2 . The result was obtained by Bardeen and Shockley 
(ref. 11). Exceptions were immediately found to the predicted T -3 / 2 variation. It was 
determined that the band structure would have to be taken into account. Herring (ref. 12) 
developed the transport properties for the ' ’many- valley” model for which the band edge 
occurs at a number of equivalent points in k- space and for which the surfaces of constant 
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TABLE I. - SCATTERING PROCESSES, VALUES OF a a , AND REFERENCES 


Types of scattering 

Constant relating 
Hall coefficient 
and carrier con- 
centration, 

Reference 

Acoustical mode (acoustical lattice vibrations of non- 
polar crystals) 

3tt/8 

11 

Piezoelectric (acoustical lattice vibrations of polar 
crystals) 

1.10 

20, 7 

Optical mode (optical lattice vibrations of nonpolar 
crystals) 


21, 22 

Polar (optical lattice vibrations of polar crystals) 

1.10 to 1.34 

23, 7 

Intervalley 


12 

Ionized impurity 

1.18 to 1.70 

24, 25 

Neutral impurity 



1.0 

i 

13 

Combination of acoustic and ionized impurity scattering 

_ . 

24, 26 

Structure defects 

1.34 to 1. 40 

1 

Electron- hole 


1 

Electron-electron 


1 


is assumed here that the effective mass part of a equals 1. 

T_ 

^Values of a are from van Daal (ref. 7). 


Direction of propagation 


t i 

O O 0-0- 

Acoustic-transverse Acoustic-longitudinal 


I 

o o -oo- 

I 

Optical-transverse Optical-longitudinal 

Figure 5. - Acoustic and optical vibrations for both 
transverse and longitudinal cases. Direction of 
wave motion is assumed to be to the right. 
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energy are multiple ellipsoids centered on each of these points. The special case of sili- 
con is illustrated in figure 4. Since the energies of conduction electrons will be within 
approximately kT of each other, those electrons which are scattered without a change in 
energy will be represented by another point on the constant- energy surfaces of approxi- 
mate thickness kT. When the scattered state is on a different ellipsoid than the original 
(such as A — C), this is called intervalley scattering. Intravalley scattering (A — B) oc- 
curs from one state to another state on the same ellipsoid. Intravalley acoustic-mode 
scattering still leads to mobility variation as T~^, but intervalley scattering can cause 
a more rapid temperature variation. Intervalley scattering is usually considered as a 
separate type and not just a variation of lattice scattering. 

The presence of impurities in the semiconductor also contributes to the scattering of 
the charge carriers. Scattering by ionized impurities yields a mobility which varies as 
T^ (refs. 24 and 25). An expression for the mobility which results from neutral impurity 
scattering was obtained by Erginsoy (ref. 13), It is not temperature dependent. 

Putley (ref. 1) outlines a method of separating the lattice scattering from impurity 
scattering. A plot of the measured mobility at a fixed temperature for a number of sam- 
ples against impurity concentrations tends to an upper limit, which indicates the lattice 
mobility. Putley states that application of this method to materials in their early states 
of development has led to incorrect results because the total impurity content is reliably 
known. 

Electron- electron, electron-hole, and structure defect scattering processes are 
briefly discussed by Putley and in detail by Blatt (ref. 6). Electron- electron scattering 
can influence other types of scattering and thus reduce the mobility. 


Suggested References 


An excellent reference is chapter 5 of reference 5. Smith gives a well-organized 
discussion of transport phenomena and scattering mechanisms. Blatt (ref. 5) discusses 
the theory of mobility for both metals and semiconductors. Putley (ref. 1) describes 
qualitatively and simply the scattering involved in nonpolar semiconductors and tabulates 
the mobility equations and gives references to their origins (pp. 146 and 147). Informa- 
tion specifically concerning silicon carbide is given by van Daal (ref. 7), and a good gen- 
eral discussion of mobility theory is presented by Ziman (ref. 8). 


DEVELOPMENT OF CHARGE CARRIER CONCENTRATION EQUATIONS 


A measurement of the Hall coefficient is a measurement of the charge carrier con- 
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centration at a specific temperature. Thus, if the carrier concentration equations can be 
developed as a function of temperature and other variables of interest, valuable informa- 
tion can be deduced from the Hall coefficient data. These equations are written for both 
the extrinsic and intrinsic cases. The notation and most of the material used herein are 
taken from chapter 3 of reference 2, Note that E v , E c , and E^ represent various en- 
ergy levels in the usual energy band scheme. However, and E a represent the en- 
ergy difference between these energy levels and the nearest band edge, that is, the con- 
duction and valence bands, respectively, and Eg is the energy difference between the 
conduction and valence band edge. 


Extrinsic Equations 


Impurity level spin degeneracy . - Consider a semiconductor with donor impurity 
atoms per unit volume which gives up one electron per atom. Each atom will be either 
ionized or neutral. Since the Fermi-Dirac distribution function f is the probability that 
any state of energy E will be occupied by an electron, the ratio of ionized donor atoms 
to neutral donor atoms is 


N 

N 


di 

dn 


1 - f 


+ E j — E 
i d c 


exp| 


kT 


(30) 


where E^ is the donor ionization energy, E^ is the Fermi level, and f is the Fermi- 
Dirac distribution function: 


f(E) = 



(31) 


Equation (30) is not quite correct. The effect of impurity level spin degeneracy must 
be taken into account by multiplying the denominator by 1//3. For the case of donor 

impurities usually /3 = 1/2, because the energy level of each atom can be occupied by an 
electron in two ways. For acceptor impurities /3 = 2, most likely. Other effects can in- 
fluence the statistical weights of the neutral and ionized conditions of the atoms. There- 
fore, /3 is used as a general symbol to include all the effects. A brief discussion of the 
effects that influence /3 is presented by Blakemore (ref. 2, p. 118). 
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Since = N di + N dn , equation (30) can be rewritten to give the density of neutral 
donors: 


N dn = 


1 + j3 exp 



(32) 


This equation is basic to the developments that follow. 

Model for single monovalent donor impurities . Figure 6 illustrates the semiconduc- 
tor model to be considered first, an n-type semiconductor with monovalent, uncompen- 
sated donors of density N^. It is impossible to have a semiconductor which would com- 
pletely satisfy this model. However, some important equations can be developed at this 
point and easily generalized later. 

The concentration of electrons excited to the conduction band is the first quantity that 
must be calculated. It is equal to the integral from E„ to infinity of the product of the 

Ks 

Fermi-Dirac distribution function and the density of electron states g(E): 


n 0 = f f(E)g(E)dE 


(33) 


More precisely, g(E) is the density of states per unit energy interval in the conduction 


* 



Conduction 

band 


N d donor level 


Figure 6. - Model of n-type semiconductor with 
monovalent impurities and no compensation. 
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band and is given by Blakemore as 


g(E) = 47 t) 



(E-E c ) 


1/2 


(34) 


where h is the Planck constant. Equation (33) is then calculated to be 

n o ■ N c r l/2<’') 


where 


N c = 2| 


( 27011 *kT 
n 


,3/2 


(35) 


(36) 


and 


r\ = 


V E c 

kT 


(37) 


and Fjygfa) is a Fermi-Dirac integral of order 1/2, Since N ( j = N ( ji +N dn and 
n Q = N^., the expression for n Q becomes 


N, 


n o = 


S -1 


1 + fT exp(e d + ?j) 


(38) 


where 


e d 



The quantity N is the effective density of states in the conduction band, and there- 
c 

fore m* which appears in N , is the density- of- states effective mass. It is not, in 
n o 

general, equal to the conductivity effective mass defined in the section MOBILITY AND 
ITS INFLUENCE. But like the conductivity effective mass, it is a function of m£, m|, 
and m|. For materials having symmetric conduction bands, the density- of- states ef- 
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fective mass is (m|m|m|)^. There is a corresponding effective density of states in the 
valance band N y and density- of- states effective mass for holes m*. 

Equation (38) simplifies if conditions are nondegenerate. If this is so, Ejygfa) and 
f(E) are approximated by expfa). A semiconductor is nondegenerate if the doping im- 
purity concentration is not too high, that is, if E f is several kT less than E c> In a de- 
generate system, only a small percentage of the electrons are capable of changing their 
energies in infinitesimal amounts. The system considered herein is nondegenerate in that 
every electron has a full opportunity to readjust to thermal change. Therefore, = 

exp ( 77 ) = n /N , and equation (38) becomes 
(J c 


2N, 


v= 


1 + 



exp(e d ) 


1/2 


(39) 


At low temperature (where kT « E d ) equation (39) can be approximated by 


E. 


n * M 1 / 2 expf-^- 
° c d V 2kT 


(40) 


Plotting In n Q as a function of 1/T results in a curve with slope E d /2k. By measuring 
this slope, E d can be approximated. Because of the dependence of N c on temperature 
(N c oT^/ 2 ) the actual slope is not exactly E d /2k. Note that equation (40) can be applied 
only if 

(1) Nondegenerate conditions apply. 

(2) The low-temperature condition prevails, that is, n Q « N^. 

(3) Little or no compensation is assumed (as for this model), that is, N a « n Q . 

A knowledge of the temperature dependence of 77 , as well as of n Q , for this model is 
valuable. Knowing 77 allows the data interpreter to determine whether the nondegenerate 
approximation can be used to simplify the carrier concentration equation. If, at a partic- 
ular temperature, 77 is less than -2 or -3, the semiconductor is nondegenerate at that 
temperature, and the simplified equation can be applied. This is the same as saying that 
the Fermi level E f is 2 kT or 3 kT below the conduction band energy. 

The variable 77 can be written for nondegenerate conditions as 



(41) 
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Substituting equation (39) into equation (41) results in 



For large T (where kT » E^), the term containing the exponential can be neglected. 
Thus, equation (42) reduces to 



(43) 


For small T (where kT « E^), the exponential dominates, and equation (42) becomes 


V = 



(44) 


It is apparent that 7 ] as a function of T starts low, rises to a maximum at 
temperature T m , and goes down again, as shown in figure 7. Therefore, if r? m is less 
than -2 or -3, the semiconductor is nondegenerate at all temperatures. There will be 
cases when the semiconductor is degenerate throughout a certain temperature range. 

This temperature range is centered around the temperature T . Thus, it is helpful to 
know not only 17 , but T m as well. Blakemore has developed a way of obtaining 
and T , and his results are summarized in the following paragraphs. 



Figure 7 . - Nominal temperature 
dependence of variable r\ for 
model of figure 6. (77 m is not 
necessarily less than r\ =0 in 
all cases.) 
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The maximum 7 } m can be calculated as a function of a quantity Z, which is itself a 
function of N^, /3, m*, and but not of T. Taking the derivative of 77 from equa- 
tion (42) with respect to T and setting the result equal to zero yield the equation 




(2e 


dm 





(45) 


where 


‘dm 


kT 


m 


z = V 


- 1 / 


in 


v3/2 


im*E d> 


(46) 


The quantity m Q is the mass of the free electron. 

A universal plot can now be made of e ( j m against Z. When this plot and a knowl- 
edge of the variables of a specific semiconductor sample which make up Z are used, the 
temperature at which 77 is maximum can be determined and then 77 m itself. This plot 
is presented in figure 8 . 

When equation (38) is differentiated and the condition dTj/dT = 0 is imposed, a 
simple relation is obtained, which when combined with the root of equation (45) leads to a 
universal graphical relation of 77 m to Z. This allows 77 m to be obtained directly from 
the curve if N^, /3, m*, and E^ of a specific semiconductor sample are known. This 
graph is figure 9. Figures 8 and 9 are reprinted with permission from F. S. Blakemore, 
Semiconductor Statistics, copyright 1962, Pergamon Press Inc. 

If there is degeneracy over part of the temperature range, equation (38) cannot be 
solved simply (i. e. , Fjygfa) * ex P( 7 ?)) - ">7 <1.3, F 1 ^ 2 ( T ?) can be approximated by 

(C + exp(- 77 )) - *, where C » 0. 27 and the electron concentration is calculated to be 


n o = 


2NjN 
d c 


< N c + cN d } + K - cN d )2 + ^‘Vc ex ^ e ^\ 


1/2 


(47) 


For 7j > 1. 3 numerical or graphical methods must be used to solve equation (38). 
Blakemore indicates the procedure to be used. 
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Model for partly compensated donor impurities . - This model assumes a donor den- 
sity compensated by a net density of acceptors N a , as shown in figure 10. This 
model also is not completely correct; several influencing effects are discussed in the fol- 
lowing sections. 

Generally, the concentration equations are similar to those for the simple model, but 
they take into consideration the presence of acceptor impurities which cause the semicon- 
ductor to be compensated. It is assumed that the donor energy levels lie at an appreciable 
level below the bottom of the conduction band and that the net acceptor levels lie even 
further away than the donor levels. The Fermi level is assumed to be well above the 
acceptor energy levels. 

The acceptors, however, acquire electrons from the donor impurities. Thus, there 
are - N a electrons to distribute between the donor levels and the conduction band. 

The number of ionized donors varies from N & at low temperatures to at high tem- 
peratures, but n Q can become larger than - N a only if the falling Fermi level ap- 
proaches another set of impurity states or its intrinsic position. 

Obviously, for this particular case, the density of conduction electrons n Q plus the 
density of electrons lost to the acceptors N & plus the neutral donors equals the 
total number of donors (i, e. , n Q + N & + N^ n = N^). Therefore, when equation (32) is 
used, 


n o +N a r- 

1 ■+ P exp (e, + v) 



Figure 10. - Mode! of compensated n-type semicon- 
ductor with monovalent donor impurities. 
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(35) 


In general, when n Q must be written as 


and A is defined as 


n o ' N c F 1/2<’>) 


' 2 ™S E d 


3/2 


A = 


h 4 


N - 


equation (48) becomes 


(49) 


-1 3/2 

1 + 0 1 exp(e d + ri) ej/ 


(v) + 


N, 

N- 


a 


(50) 


# 


When ri < 1.3, 
duces to 


^1/2 ( ? ?) can approximated by (C + exp(-?])) - *, and equation (50) re- 


2N c (N d -N a ) 


[n c + C(N d - N a ) + /3 _1 N a exp(e d )] + -fa, - C(N d - N & ) + /3 _1 N a exp(e d )] 2 + 4^ -1 (N c + CN a )(N d - N a )exp(e d ; 


1/2 (51) 


When the Fermi level is several kT below the conduction band (i. e. , t] = -2 or -3), 
F l/2^) “ ex P(^)» and equation (50) reduces to 


n = 
o 


2 ( N d - N a ) 


N 

1 +— exp(e d ) 

0N„ d 


r r 


+ < 


l +■ 


i3N, 


exp(e d ) 


VT 


1/2 


(52) 


i3N, 


( N d - N a )exp(e d ) 


Both equations (51) and (52), when N a goes to zero, reduce to equations (47) and 

(39), respectively, which were originally derived for the uncompensated model. 

When compensation is small, there exists a temperature range where 

N « n « N .. Equation (52) becomes, in this case, 
a o ci 
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(53) 


n o = (^ N c N d )1/2 ex P 



which is the same as equation (40). At lower temperature, however, when n < N , 

O a 

equation (52) becomes 


<54> 

Note the dependence of In n Q on 1/T. The dependence is twice as large for equation (54) 
as for equation (53). Since In n Q as a function of 1/T is often plotted to determine 
from the slope of the curve, care must be taken to differentiate between these equations. 
When compensation is large, equation (53) is never valid. For heavy compensation, 
equation (52) becomes, at all temperatures, 


n o 


1 +■ 


N,, - N q 
d a 

N„ 


/3N, 


exp(e d ) 


(55) 


Obviously, the equations for single monovalent and partly compensated donor impurities 
must be applied cautiously to a real semiconductor. Forethought is needed to prevent 
serious errors when the experimental data are interpreted. 

Model for compensated semiconductor with two monovalent donor impurity levels . - 
This model assumes two monovalent donor impurity levels which are compensated by a 
net acceptor density. The model is illustrated in figure 11. The concentrations of the 
two donor levels are denoted as and respectively, and their energies are 

written as and E^ with E^ > E^g* 

The acceptor impurities trap N„ electrons so that the net density available is 

a 

N dl + ^d2 “ ^a" T ^ e r Native magnitudes of N^, N^, and N a must be considered, 
while remembering that E^j > E^ is assumed. It is found that only two types of be- 
havior result: N a > and N a < N dr K N a > N^, there will be no electrons at the 
Ndi level at all temperatures. This corresponds to a semiconductor which has a single 
type of donor impurity with an effective compensating density of N a - N^. Therefore, 
the analysis for partly compensated donor impurities is applicable. 

When N a < N^, the levels will be full, and N^j - N a electrons will be in the 
N dl * eve * states at 0 As the temperature rises, the N^ donors start giving their 
electrons to the conduction band. When all the N^ donor levels are empty, the elec- 
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Conduction 

band 


% donor level 
Nd 2 donor level 


N a net acceptor 
level 


Figure 11. ~ Model of compensated n-type semiconductor 
with two independent monovalent donor impurities. 



Figure 12. ~ Action of electron concentration as function of temperature for compensated 
semiconductor with two independent monovalent donor impurities. 


trcm concentration levels off. As the temperature increases it begins to rise until all the 
N d 2 electrons are ionized, and then it levels off again. The electron concentration as a 
function of temperature is shown in figure 12. 

it. 

To discuss the model mathematically, is defined at the j ■ monovalent donor 
impurity (where there are a total of M different types) with ground state binding energy 
E^. Since the density of ionized donors is equal to the conduction band electron concen- 
tration plus the density of acceptors, the following equation holds: 


28 



M 



N dj 

1 + 07 1 exp(r? + e d -) 


( 56 ) 


And if the semiconductor is nondegenerate, exp ( 77 ) can be approximated by n /N . 
Equation (56) can then be analyzed. 

There are other extrinsic doping models that are important but will not be presented. 
Some of these include multivalent doping atoms, amphoteric impurities (i. e. , an impurity 
that can act either as an acceptor or a donor), and those cases influenced by lattice de- 
fects. For more information about these types of semiconductors see Blakemore (ref. 2). 

Influence of excited states on carrier concentration equation . - An analogy exists be- 
tween the energy level structure of atoms and the energies that a donor impurity may ex- 
hibit. Not all donor electrons have an energy at the level E c - E^. Rather this level is 
the donor electron ground state, and the donor electron could be in an excited state at an 
energy above the donor ground state. Figure 13 depicts an extended n-type compensated 
semiconductor problem, which is considered in the following paragraph. 

Let r denote the energy state. For example, r = 1 is the donor ground state, and 
r = 2 is the first excited state. Then each state has a spin degeneracy /3~* and lies at 
energy E^ = kTe^ above the ground state. The donor can bind an electron in a state 
only if it does not also have an electron bound in another state. 



Figure 13. - Compensated n-type semiconductor of figure 10 
but including donor excited states. 
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The resulting equation for the nondegenerate case is 


n o = - 


N 

1 + — — exp(e d )(l + S) 

N A 


2(N d - N a ) 


N 

1 + — exp(€ d )(l + S) 

N A 


2 1/2 
2 4(Nj - N ) 

+ — exp(e d )(l + S) ' 

Vl ' 


( 57 ) 


where 


s ^ZS exp< ' si) 

r=2 


(58) 


Equation (57) is quite similar to equation (52). If S = 0, they would be identical. Thus, 
the equations previously derived are influenced by the action of excited states in a simple 
manner. To proceed from one equation to another the term 1 + S must be appropriately 
substituted. Note that at low enough temperatures (where kT « E^), S « 1 and 
equations (57) and (52) coincide. At high temperatures, n Q = N R . Thus, only at interme- 
diate temperatures does S affect n Q significantly. 

Other factors which influence carrier concentration equations . - It has been found 
that the ground-level donor energy is sometimes split into two levels which are A apart. 
This effect is a function of energy band structure. Long and Myers (ref. 14) have meas- 
ured A for a specific case. The splitting influences the electron or hole concentration, 
but the details and the degree of influence are not discussed herein. 

Strain (ref. 15) also affects the structure of impurity states in most materials as 
well as a magnetic field (ref. 16). However, except for large strains or large magnetic 
fields, these effects are negligible. 


Intrinsic Equations 

A semiconductor displays intrinsic conduction when the temperature is high enough 
that the carrier concentration due to thermal excitation from the valence band is larger 
than the concentration due to excited impurities. Analysis of data in the intrinsic region 
yields the band gap of the semiconductor. 

Before the carrier concentration equation is written, clarification is needed concern- 
ing the total hole density p Q . The equation for total electron density n Q has already 
been written. The hole density can be expressed similarly as 
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p 0 te Vi/a<-V>> 


(59) 


where 


N y = 2| 


'2?rm*kT 


,3/2 


(60) 


and £g is the band gap in units of kT(€ g = E^/kT). The hole density p Q is the integral 
over energy from minus infinity to the top of the valence band E y of the product of the 
density of states and the Fermi distribution for holes, which is 1 minus the Fermi distri- 
bution for electrons. 

In most cases (especially wide-band-gap material) the nondegenerate approximation 
is valid because the gap is a fairly large multiple of kT and the Fermi level is near the 
center. Since in an intrinsic semiconductor n. = n Q = p Q (where m is defined as the 
intrinsic concentration of holes and electrons), 


n i = N c exp ^ = N v ex P(" e g _1? ) (61) 

The intrinsic carrier concentration equation is obtained when the extreme left side of 
equation (61) is multiplied by n^, the extreme right side is multiplied by N c exp fyj), and 
the square root is taken. 


n. = (N N J 1 / 2 exp(-i' 
1 v c \ 2 kT ; 


\ 


Another basic equation can be derived from equation (61): 


Vo = n i 


(62) 


(63) 


INTERPRETATION OF DATA 

Measured values of Hall coefficient and resistivity against temperature are usually 
displayed in a standard way. The carrier concentration n Q is calculated from the Hall 
coefficient, and then the natural logarithms of n Q and p are plotted as a function of in- 
verse temperature, as in figure 14. (An n-type semiconductor is assumed throughout 
this section. ) Each part of the curve on each plot yields a particular bit of information 
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Figure 14. - Typical plots of natural logarithms of carrier concentration and of resis- 
tivity as function of inverse temperature. These plots only indicate general shape 
of curves. 


about the semiconductor. The general shapes of the curves are shown in figure 14. Note 
that the slopes of the curves at both extremes of temperature are proportional to either 
or Eg. In the extrinsic region, the carrier concentration rises as the temperature 
increases because more and more electrons are being excited into the conduction band. 
When all of the available electrons from the doping atoms are ionized, n Q will reach a 
constant value (N^ - N a if compensated). As the temperature rises, hole-electron pairs 
are generated as the thermal energy available becomes comparable to E g . Thus, n Q in- 
creases rapidly. 


32 



Equations have been developed for each part of each curve in figure 14. The inter- 
pretation of the data involves the choosing of the correct parameters in the equations. 

The correct parameters, when substituted in the appropriate carrier concentration equa- 
tion, will match the calculated results to the experimental data. The correct parameters 
are obtained by an iterative process. 


Extrinsic Data 


Several extrinsic equations have been developed in the previous sections for degener- 
ate and nondegenerate semiconductors. Some estimate of the degeneracy must be made 
in order to pick the applicable equation. The degeneracy is determined by the position of 
the Fermi level. Since the Fermi level varies as a function of temperature (as well as 
doping), it is possible for a particular equation to be good in one temperature range and 

not good in another. Normally, a semiconductor must be doped to 10*® or 10^ centi- 
-3 

meter before it becomes degenerate. Therefore, the procedure is to estimate the max- 
imum value of the Fermi level. And if it is several kT below the conduction band, the 
semiconductor is always nondegenerate. If it is at or above the conduction band, the 
Fermi level must be calculated more carefully as a function of temperature. 

Blakemore' s variable 77 as a function of temperature yields the degeneracy informa- 
tion. Even for a compensated semiconductor, 77 m and T m (where T m is the tempera- 
ture at which 77 = 77 m ) can be determined reasonably well by substituting N d - N a into Z 
(eq. (46)) and then reading T m and r} m from figures 8 and 9. If equation (52) is found 
to be unusable, equation (50) or (51) must be utilized. 

Hall extrinsic data . - For this discussion it is assumed that the semiconductors in 
question are nondegenerate. However, much of the material is still applicable to the de- 
generate cases. 

Blakemore (ref. 2) and Putley (ref. 1) both outline a method of interpreting experi- 
mental Hall coefficient data which Blakemore calls one of the most popular approaches. 
The method starts with a rearranged form of equation (52) 


n o< n o 4 N a> 
(N d - N a - n 0 ) 




(64) 


where the model for partly compensated donor impurities is used. The effects of excited 
impurity levels are neglected. Notice that the right-hand side depends on the character 
of the impurities and on the temperature, but not on the acceptor or donor densities. The 
first problem is to find values of N a , N d , and E^ which satisfy the equation at all tem- 
peratures. Each of these parameters can be estimated. Then, by iteration the param- 
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eters are adjusted until the calculated carrier concentration that results (using N a , N^, 
and E^) matches as closely as required with the experimental carrier concentration data. 
Blakemore states that the method of least squares applied to equation (64) will give those 
values of N a and which give the best fit over the entire temperature range. Satis- 
factory results are obtained by using a repetitive method of varying the parameters until 
the calculated results match the experimental data (private communication from Robert P. 
Ulman, formerly of Lewis Research Center and now at the Addressograph-Multigraph 
Corp. of Cleveland, Ohio). These curve-fitting procedures are time consuming, and dig- 
ital computers are helpful in doing these calculations. 

There are several ways of obtaining the estimates of the parameters in question (as- 
suming an n-type semiconductor). Blither /3 or m* (which is in the equation for N ) 
must be known independently. That is, if m* is known by some other type of measure- 
ment, such as cyclotron resonance, N c can be calculated. The quantity - N a can be 
estimated from the plateau (exhaustion region) of the curve of In n Q against reciprocal 
temperature. Donor ionization energy is obtained from the slope of the In n Q 
against 1/T curve at low temperature (see eqs. (53) and (54)). Finally, N a /j3 can be 
calculated by using equation (54) if n Q at an arbitrary low temperature is known. 

Putley suggests two methods for estimating and N a . In the first method the in- 
tercepts of the low-temperature part of the plot of In n Q against 1/T, where equa- 
tion (54) is valid, are measured. This allows (N^ - N a )/N a to be determined. Then, 
from the exhaustion region of the same plot, - N a is estimated and, therefore, values 
of and N a are obtained. In the second method, the curve at low temperature in fig- 
ure 15 is extrapolated backwards using equation (54) until it intercepts the exhaustion line. 
At the temperature corresponding to the point of intersection, the left-hand side of equa- 
tion (64) equals N a . By estimating - N a from the exhaustion region, and N a 
can be determined. 

Often there will occur plots of In n Q against 1/T in which there is more than one 
dip in the curve. This situation might correspond to semiconductors with two or more 
levels of either donor or acceptor impurities, or to impurities which can accommodate 
more than one electron, or to structural defects of the lattice of the semiconductor. 

When this is the case, the generalized equation (56) must be used in the interpretation of 
the data. Putley includes an excellent discussion of this problem (ref. 1, section 4. 8). 

Because of the temperature dependence of the quantities N c and N y , the plot of 
In n Q against 1/T will not be a perfectly straight line at low temperature. This makes 
the calculation of the impurity ionization energy E^ inaccurate. However, if the tem- 
perature dependence of N is plotted on the ordinate along with n , the slope at low 

c o 

temperatures will be directly proportional to E^. For example, assume that conditions 
are such that equation (54) holds at low temperatures. Notice that N c is proportional 
to Therefore, ln(n Q /T^) can be plotted against 1/T to obtain a curve with a 
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Figure 15. - Method of determining concentration of donor and acceptor impurities and 
N a of compensated semiconductor. 


slope truly proportional to E^. 

It must be emphasized that the interpretation of the data need not follow the pattern 
outlined herein. There are other variations which have been used and published. For ex- 
ample, van Daal, Knippenberg, and Wasscher (ref. 17) have developed a method for anal- 
yzing Hall data which extrapolates the measured curves to the region of complete exhaus- 
tion. This method could be valuable in those cases where data could not be taken at suf- 
ficiently high temperatures. 

The fact that impurity ionization energy becomes smaller as the impurity density is 
increased was shown by Pearson and Bardeen (ref. 18). At a certain density, character- 
istic of the host lattice and the impurity type, the impurity ionization energy actually be- 
comes zero. The variation in E„ as a function of N is shown in figure 16. Care must 

a a 

be taken when determining the impurity ionization energy when the doping is high. 

Resistivity extrinsic data. - Little has been said about the interpretation of resistiv- 
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Figure 16. - Effects of boron concentration on boron acceptor ionization energy in silicon 
according to measurements of Pearson and Bardeen (ref. 18). 


ity against temperature data because only the impurity ionization energy can be deter- 
mined from it. However, resistivity is simpler to measure than the Hall coefficient. 
Therefore, in the case where the Hall data (In n Q against 1/T curve) does not straighten 
out at low temperatures and, thus, the ionization energy is still in doubt, resistivity data 
may be very useful. In this case, since 1/p = n o qp n (for n-type) and the mobility is not 
an exponential function of T at low temperatures, the exponential exp(-E^/kT) in the 
carrier concentration term n Q dominates the behavior. A plot of In p against 1/T 
yields the ionization level. Resistivity measurements are then emphasized at lower tem- 
peratures, and a more accurate determination of the impurity ionization energy may be 
possible. 

Resistivity, as well as the Hall coefficient, is required to calculate the Hall mobility. 
A knowledge of mobility is important for device design. 


Intrinsic Data 

When a semiconductor is intrinsic, n Q = p Q = n., and the equations of resistivity and 
Hall coefficient when both types of carriers are present reduce to 
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and 


J ‘ * a lK + V 


(65) 


iu. jL^lL3 (66) 

« n l <*p + “n 

The equations are altered from the extrinsic case by the inclusion of sums and differences 
of the electron and hole mobilities. For example, the measured Hall coefficient cannot 
be directly related to the carrier concentration, as in the extrinsic case. Thus, the car- 
rier concentration cannot be calculated unless the mobilities are known. However, the 
band gap energy Eg appears in the exponential term of equation (62) for the intrinsic 
carrier concentration. This term dominates the temperature dependence of the mobili- 
ties, and thus the slopes of the plots of In Rjj against 1/T or In p against 1/T 
(fig. 17) give the band gap energy. Once again, as in the extrinsic case, it may be easier 
to measure resistivity at the temperature extreme (i, e. , high temperature) of intrinsic 
conduction. The band gap energy Eg usually is a function of temperature (ref. 2). But if 
it is assumed to be a linear function of temperature, 

E = E - aT (67) 

g go 

the exponential term in the carrier concentration equation becomes a constant times 

exp(-E /2kT). Therefore, the plots in figure 17 yield the band gap at 0 K. The band 
go 

is not necessarily a linear function of temperature, therefore the plots still appear 
;urved in these cases. 

M the Hall coefficient is divided by the resistivity, the Hall mobility is the result, 
intrinsic conduction it becomes 


= a (M p “ M n ) (68) 

Interpretation of Data Near Intrinsic Conduction 

The Hall and resistivity curves as a function of temperature take on unexpected 
shapes as they approach the intrinsic region. These shapes are dependent on whether the 
semiconductor is n- or p-type and on the semiconductor mobility ratio b = P n /M p* ® is 
possible that the intrinsic data will be misinterpreted if these shapes are ignored. The 
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(c) For b>3. 7. (d) For all b. 

Figure 17. - Variation of natural logarithms of Hall coefficient and resistivity with inverse temperature for three mobility 
ratios, b < 1. 46, 1. 46 < b < 3. 7, and b > 3. 7. Dashed li nes i ndicate true i ntri nsic slopes. 

general shapes of the curves are shown in figure 17. Dunlap (ref. 19) first discussed 
this problem, but Hunter (ref. 9) describes it in more detail. 

There are two prominent features of the curves in figure 17. The first is that the 
Rjj for p-type semiconductors usually changes sign as the intrinsic region is approached. 
This feature is expressed in equation (6), since M n is usually greater than u . The 
second is that Rjj and p for some cases at some temperatures are greater tnan their 
intrinsic values (shown by the dashed line). This is known as overshoot. Thus, the re- 
sistivity for p-type semiconductors for all mobility ratios is greater than the intrinsic 
resistivity in a certain temperature range. As the temperature increases, the resistivity 
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approaches its intrinsic value. It is now obvious that if measurements of Rg and p are 
not made at temperatures high enough to go beyond the overshoot region, the measured 
slope will be incorrect. And, therefore, the band gap energy will be incorrect. 


SPECIFYING TEMPERATURE RANGE OF DATA 

The Hall coefficient and resistivity must be measured throughout certain specified 
temperature ranges. The informative parts of the curves of n Q and p as a function of 
1/T appear only at characteristic temperatures. Therefore, in order to ensure that 
enough data are collected, the required temperature range must be estimated beforehand. 
This is also likely to be important in terms of cost. For example, if the electrical prop- 
erties of a specimen need not be measured below liquid- nitrogen temperatures (77 K), the 
large cost and complications of a liquid-helium cryostat can be avoided. 

At low' temperatures the limiting factor is the ionization energy of the impurities. 

The plot of In n Q against 1/T must have a constant slope at low temperatures. As the 
energy levels of the impurities lie closer to the conduction band, the lowest temperature 
at which data must be taken decreases. 

The high-temperature limit is dependent on what is to be obtained. If the band gap of 
the basic material is to be measured thermally, the data must extend well into the intrin- 
sic region. To measure the doping concentrations of a particular piece of semiconductor, 
the measurement of the carrier concentration (i. e. , Hall coefficient) must extend to the 
exhaustion region. There all the impurities are ionized and the carrier concentration is 
constant with temperature. 

Obviously, the low- and high-temperature limits are peculiar to each semiconductor 
material. Therefore, the method of specifying the temperature range will involve search- 
ing the literature for information about the band gap and ionization energies and compar- 
ing this with known semiconductors, such as silicon and germanium. Measurements on 
these semiconductors must be as low as liquid-helium temperatures. 

In the case of a wide-band-gap semiconductor, such as silicon carbide, the tempera- 
ture range may be shifted upwards and expanded. Since the band gap is approximately 
3 electronvolts for hexagonal silicon carbide, the plateau region of carrier concentration 
data will probably be above room temperature. Previous measurements made at the 
Westinghouse Research Laboratory in Pittsburgh indicate that temperatures of at least 
550 to 600 K are required to obtain the necessary extrinsic data. This is nearing the in- 
trinsic region, which occurs at temperatures above 850 K (for doping levels of about 
10 17 cm" 3 ). 

The low -temperature limit for silicon carbide can be as low as liquid-helium tem- 
peratures, but, practically speaking, liquid- nitrogen temperatures are sufficient. An 
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impurity energy level may be very close to a band edge and thus become ionized at ex- 
tremely low temperatures. However, the probability of this happening in silicon carbide 
is probably less than in silicon (which has a band gap about one-third that of silicon car- 
bide). Therefore, for the example of silicon carbide, the temperature range of interest 
would be from liquid-nitrogen temperatures to greater than 550 K. Probably the upper 

limit should be extended to 850 K to make certain the exhaustion region is reached. For 

17 -3 

doping levels greater than about 10 centimeter , the upper limit may need to be in- 
creased even further. 


CONCLUSIONS 

The measurement and interpretation of the Hall coefficient and resistivity of semi- 
conductors were outlined, and the basic equations of carrier concentration, obtained from 
other sources, were summarized. Some known methods of interpreting the experimental 
data by use of these basic carrier concentration equations were then discussed. 

In considering the applicability of the measurements to both basic materials research 
and device design, it was pointed out that the energy levels of the impurities, the band 
gap energy, and the mobility are valuable to the materials researcher. For the device 
designer, the mobility, the impurity ionization energies, and the doping concentrations 
are important. All of these quantities are obtained from an analysis of the Hall and re- 
sistivity data. 

Two methods of measuring the Hall coefficient and resistivity were briefly described. 
The van der Pauw method is favored because in that method the semiconductor sample 
need not be of a particular shape. This is advantageous in the case where only small 
samples are available. 

The mobility of semiconductors was shown to be related to an energy- averaged re- 
laxation time which is a characteristic of a particular scattering process and a specific 
semiconductor material. The results of the derivation of the relation between mobility 
and relaxation time were given for the simplest case, where the constant- energy surfaces 

j 

were spherical and there was only one energy minimum or maximum. The results were 
extended to the cases where the constant- energy surfaces were nonspherical and there 
were multiple energy minimum s or maximums. Although no relaxation time calculations 
for specific scattering processes were given, the different kinds of scattering were qual- 
itatively discussed. The factor a which relates the Hall coefficient to the carrier 
concentration resulted from the same type of derivation. 

The equations for charge carrier concentration were developed for both extrinsic and 
intrinsic conduction. The extrinsic equations were found to be primarily functions of 
temperature, impurity concentrations, and the impurity ionization energy. The intrinsic 
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equation was a function of temperature and the band gap energy. The extrinsic equations 
were found to be modified by different types of doping, by inclusion of excited states of 
the doping impurities, and by degeneracy. 

This set of equations for extrinsic and intrinsic conduction was then used as the basis 
for developing methods of interpreting the measured experimental data. It was found that 
the parameters of the equations could be estimated with varying degrees of accuracy from 
the plot of In n Q against 1/T, where n Q is the equilibrium concentration of electrons 
and T is temperature. By a method of iterative calculations, the estimated parameters 
were varied until the calculated results fit the experimental data curve. It was shown that 
the resistivity plot may help to obtain more accurate values of the impurity ionization en- 
ergy and the band gap energy at extremely low and high temperatures, respectively. 

The shape of the curves of the Hall coefficient and resistivity against reciprocal tem- 
perature near the point where intrinsic conduction begins was illustrated. It was found 
that knowledge of these curves was necessary for correct interpretation of intrinsic data. 

The temperature range over which the Hall coefficient and resistivity must be meas- 
ured in order to use the methods presented in this report was discussed. Measurements 
must be collected at a temperature low enough that the slope of the In n Q against 1/T, 
curve becomes constant. At high temperatures the same curve must level off at a con- 
stant n Q (the exhaust region). If intrinsic properties of the semiconductor, such as band 
gap energy, are wanted, the measurement must extend well above the temperature where 
the semiconductor becomes intrinsic. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 27, 1968, 

120-27-01-25-22. 



APPENDIX - SYMBOLS 


A 




E, 


E 


g 

go 



defined function not dependent 
on temperature 

magnetic field, tesla 

conductivity mobility ratio, 

effective mass ratio, m£/m|, 
energy, eV 

acceptor ionization energy, eV 

energy of bottom of conduction 
band, eV 

donor ionization energy, eV 

donor ionization energy of 
donor of j l kind, eV 

Fermi level, eV 

band gap energy, eV 

Eg at T = 0 for linear varia- 
tion of E , eV 

B 

energy difference between 

th 

ground level and r level, 
eV 


g(E) density of states per unit en- 

ergy interval in conduction 
band, cm"^ eV~* 

-14 

h Planck constant, 6.63x10 

J sec 


I 

J , J , J 
x’y’z 


current, A 

components of current density, 

A — 2 

A cm 


k Boltzmann constant 

1.38X10" 23 J K" 3 


k wave vector '(components 

ky, k z makeup 

k- space), m 


m 

mj 

m2 


m* 

n 





m. 


(effective) length of bridge- 
shaped crystal sample, cm 

components of effective mass 
tensor, kg 

effective mass of electrons, kg 

mass of free electrons, 

9. llxlO" 31 kg 


E y energy of top of valence band, 

eV 


***?** ^“Ponents ol electric field, 
Vcm 

F1/2O7) particular Fermi-Dirac in- 
tegral 


f(E) Fermi-Dirac distribution 

function 



correction factor in van der 
Pauw's resistivity equation 


m* 



N, 




effective mass of holes, kg 

concentration of acceptor im- 
purities, cm -3 

effective density of states in 

-3 

conduction band, cm 

concentration of donor im- 
purities, cm -3 

concentration of ionized donor 
impurities, cm" 3 

concentration of donor im- 

xi* ^ o 

purities of j kind, cm" 
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N 


dn 


N, 


n. 


n_ 


R 


ab, cd 


% 

S 



t 



concentration of neutral 

o 

donor impurities, cm 

W 

width of bridge- shaped crys- 
tal sample, cm 

effective density of states in 
-3 

valence band, cm 
intrinsic concentration of 

Z 

defined function not dependent 
on temperature, cm’* 

eV -3/2 

electrons and holes, cm 

equilibrium concentration of 
electrons, cm”'* 

a 

constant relating Hall coef- 
ficient and carrier concen- 
tration 

equilibrium concentration of 

_3 

holes, cm 


impurity level spin degener- 
acy 

electronic charge, 
1.60X10” 19 C 

h 

impurity level spin degener- 
acy of impurity of j** 1 kind 

ratio of voltage measured 
across contacts d and c 


impurity level spin degener- 
acy of r n level 

with current flowing into 

€ 

defined variable, (E - E c )/kT 

a and out b (other such 
ratios defined similarly), 

e d 

defined variable, E^/kT 

ohms 

e dj 

defined variable, E^. /kT 

Hall coefficient, cm* C”* 

e dm 

defined variable, /kT m 

summation of effects of en- 

e 

g 

defined variable, E /kT 

O 

ergy of excited state of 

e rl 

defined variable, E r ^ /kT 

impurity atoms 

V 

defined variable, (E^ - E c )/kT 

temperature, K 

77 m 

maximum of 7] 

temperature where r\ - q , 
K 

M 

conductivity mobility, 

2 -1 -1 
cm V sec 

thickness of crystal sample, 
cm 

m h 

m L ,m T’ \ 
M n 

Hall mobility, cm* v”* sec”* 
components of conductivity 

voltage drop from C to E, 
V 

Hall voltage, V 

2 -1 -1 
mobility, cm V sec 

electron mobility, 

2 -1 -1 
cm V sec 

x component of average 
drift velocity of electron, 
cm sec”* 


2 -t -1 

hole mobility, cm V sec 
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p resistivity, fi/cm 

a electrical conductivity, cm 


r relaxation time, sec 

(r) relaxation time averaged over en- 

ergy, sec 
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